function d1 = f1(t,x,y,z)
m = 0.43;
f_x = F_d(sqrt(x^2+y^2+z^2))/sqrt(x^2+y^2+z^2)*x;
d1 = -f_x/m;
